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Field of the invention 

The present invention relates broadly to a method of improving a grating design function 
describing a refractive index variation defining a multi*channel grating stmcture in a waveguide 
5 material, to a method of fabricating a multi-channel grating structure, and to a multi-channel 
grating structure. 

Background of the invention 

Multi-channel grating structures are typically writtm into photosensitive waveguides. 
The grating stmcture comprises a refractive index profile induced in the photosensitive 
10 waveguide, which in turn detemiines the optical characteristics such as the reflection, 
transmission, and groiq> delay characteristics of the resulting grating stmcture. 

The amplitude and phase of the refractive index profile can be described by a grating 
design function, which in tum is used in the control of an apparatus for writing the grating 
stmcture into the photosensitive waveguide. Typically, die writing apparatus comprises an 
1 S optical interferometer. 

It is known that a maximum refi^tive index contrast required to write multi-channel 
gratings increases with increasing nimiber N of channels. Since photosensitive materials, such 
as a photosensitive fibre used to fabricate a Bragg grating, have material limits of the maximum 
achievable photo-induced refractive index change A??, this represents a limitation on the 
20 maximum number of channels N that can be recorded in a given material. Thus, it is highly 
desirable to reduce a required Ans as much as possible. 

At least preferred embodiments of the present invention seek to provide a method of 
improving a grating design function in terms of a reduced maximum refiractive index variation 
required in the waveguide material along the grating stmcture, without compromising the 
25 quality of spectral characteristics. 

Summary of the invention 

In accordance with a first aspect of the present invention there is provided a method of 
improving a grating design function describing a refractive index variation defining a multi- 
channel grating stmcture in a waveguide material, the improvement being a reduced maximum 
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refractive index variation in ttie waveguide material along the grating structure while 
maintaining a desired functional spectral domain in a spectral response function associated with 
the design function, the method comprising the steps of: 

- modifying a first design function to generate a second design function having a reduced 
5 maximum amplitude compared witii the first design fimctioxi, 

- determining a second response function associated with the second design function, 

- modifying the second response function to create a third response function having a 
desired functional spectral domain, and 

• determining a third design function associated with the third response function, 

10 and iterating the method steps until the desired improvement is achieved, wherein the 

thud design function of the previous iteration takes the place of the first design function of the 
next. 

In one embodiment, the step of modifying the second response function comprises 
replacing the corresponding spectral domain of the second response function by the desired 
1 5 functional spectral domain. 

The desired functional spectral domain may conq)rise a corresponding spectral domain 
of a first response function associated with the first design function. 

The method may comprise the pre-iteration step of determining the first design function 
fi-om a or the associated first response function. 

20 Preferably, the det^mining of the response functions fix>m tibte design functions 

comprises solving a direct scattering problem, and tiie determining of tiie design functions tcom 
the response functions comprises solving an inverse scattering problem. 

In one embodiment, the step of modifying the second response function to create the 
third response function is conducted in a maimer such that the third response function has a 
25 desired response characteristic in at least one spectral domain other than the functional spectral 
domain. The step of modifying the second response function may comprise replacing the 
corresponding at least one spectral domain with the desired response characteristic. 

Preferably, the step of modifying the first design function comprises applying a 
normahsation process or an averaging process to the first design function. 
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In one embodiment, the noimalising process comprises replacing an amplitude function 
k{z) of the design function by a product of the square root of a constant A and a corresponding 
single-channel seeding ampUtude function k^{z), while maintaining a phase function of the 
design function. The constant A may be deiSned by a normalisation condition. In one 

5 embodimait, A is defined by the normalisation condition ^K^dzl ^K]dz^ where / is a 

0 0 

length of the multi-channel grating structure. 

In one embodiment, the averaging process comprises averagmg over a san^ling or 
quasi-sampling period of the design fimctioa 

In another embodiment, the normalising or averaging processes n^iay be complimented or 
10 replaced by a nonlinear transform reshaping ("squeezing") operation, which reduces k^{z) 

while keeping the parameter Jir^db substantially unchanged. 

0 

In accordance with a second aspect of the present invention, there is provided a method 
of fabricating a multi-channel grating structure, the method comprising the step of improving a 
grating design function describing a refractive index variation defining the multi-channel 
15 grating stmcture in a waveguide material as defined in the first aspect. 

In accordance with a third aspect of the present invention, there is provided a multi- 
channel grating structure £abricated utilising a method of fabrication as defined in the second 
aspect 

Brief description of the drawings 

20 Preferred forms of the present invention will now be described with reference to the 

accompanying drawings. 

Figures 1 (aHd) show the calculated first spectral characteristics and the first design 
ampUtude of a multi-channel grating structure (the first design phase is not shown), illustrating 
one stqi of a method embodying the present invention. 

25 Figures 2 (aH^) show the calculated second spectral characteristics and the second 

design results (normalised/averaged amplitude profile), illustrating one step of a method 
embodying the present invention. 
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Figures 3 (a)-(d) show the calculated third spectral characteristics and the third design 
results which correspond to the first completed iteration embodying the present invention. 

Figures 4 (aHd) show calculated spectral characteristics and design results after 8 
iterations of an improvraient process embodying the present invention. 

5 Figures 5 (a)-(d) show calculated spectral characteristics and design results after 150 

iterations of an improvement process embodying the present invention. 

Figure 6 is a flowchart illustrating a multi-channel grating desi^ improvemmt process 
^bodying the present invention. 

Figure 7 shows an experimental set-itp for writing a multi-channel grating structure of a 
1 0 multi-channel grating design embodying the present invention. 

Detailed description of the embodiments 

The preferred onbodiment described provides a method of in^roving a grating design 
function in temis of a reduced maximum refiractive index variation required in a waveguide 
material along the grating stracture. 

15 Figures l(a)-(d) show calculated spectral-characteristics and design of a 17-channel 

Bragg grating prior to iteration procedure. It will be ^predated by a person skilled in the art 
that the spectral-charactmstics are of a supeifo quality with virtually no deviations firom a 
desired spectral-characteristics, in the example shown in Figure 1 a square-sh^ed transmission 
characteristic (compare Figure lb). 

20 It has been found by the inventors that the design of the 17-channel Bragg grating shown 

in Figure 1 can be inq)roved in terms of minimisation of a maximum refractive index change 
required in the photosensitive material, e.g. a photosensitive optical fibre, into which the 17- 
channel Bragg grating is to be written. In other words, the fibre Bragg grating (FBG) amplitude 
versus distance z along the FBG (Figure la), can be improved in terms of mmimisation of the 

25 maximum amplitude without compromising the spectral response. In the example embodiment, 
this is achieved by subjecting tiie design to an iterative optimisation procedure which will now 
be described. 

In the example embodiment, the maximtnn in the FBG ampUtude versiis distance z is 
initially reduced in the following way: Starting from the FBG amplitude versus distance z 
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shown in Figure 1(a), the amplitude function k{z) is replaced by V^c, (z) , where k^{z) is the 
coiresponding single channel (seeding) grating an:5)litude function and Ae constant A is defined 

by the normalisation condition A = ^K^dzl ^K]dz , Alternatively k{z) is replaced by its 

0 0 

averaged (smoothen) version. Alternatively, the normalising or averaging process may be 
5 complimented or replaced by any nonlinear transform resha5)ing ("squeezing") operation, which 

reduces (z) while keq)ing the parameter ^K^dz unchanged. The unchanged mxilti-channel 

0 

grating phase and the modified (the second) grating ampUtude are then used as the input data for 
solving a direct scattering problOTi utilising a suitable direct scattmng solver algorithm. Thus, 
a spectral response function corresponding to the modified grating design function is 
10 determined. 

The second grating design amplitude and the coiresponding spectral response data are 
shown in Figures 2(a-d). Clearly the spectral characteristics in the central part 100 of the 
wavelength range m Figure 2(b), which is the fimctional spectral domain in the example 
embodiment, are less than perfect (compare central part 102 in Figure 1(b)). As a next step these 
IS characteristics are modified. 

Within flie central channel spectral domain it is replaced by the initial (pre-iteration) 
data, compare central part 102 in Figures 1(b) and 3(b). Outside tixe cratral channel spectral 
domains the determined spectral response data is left intact, compare regions 104 and 106 in 
Figures 2(b) and 3(b), i.e. that spectral domain is allowed to evolve during the iteration process. 

20 As a final step of each iteration, an mverse scattering problem is solved for the modified 

spectral response data to determine the associated, better optimised, multi-channelled gratmg 
design function, see Figures 3(a) and (c). 

The above method steps are then iterated in the example embodiment, until a desired 
improvement is achieved. Figures 4 and 5 show the level of convergence after 8 and 150 
25 it^tions for the 17-channel grating design. For each of the Figures 4 and 5, (a to d) show 
calculated spectral characteristics and design of the 17-channel FBG. As can be seen from a 
comparison of Figure 1(a) and Figure 5(a) the optimisation process embodying the present 
invention has improved by about 13.5% to a lower value, while, due to the nature of the 
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aptimisation process embodying llie present invention, the superb spectral characteristics of the 
design (compare Figures 1(b), (d), and Figures 5(b), (d)). 

It will be ^preciated by a person skilled in the art that the step of modifying the second 
response function to create tiie third response function may further be conducted in a manner 
5 such that the modilSed response function has a desired response characteristic in at least one 
spectral domain other than the functional spectral domain, by e.g. replacing the corresponding 
at least one spectral domain in the second response function with the desired response 
characteristic. At Hxe same time, it will be appreciated that parts of the spectral domain need to 
remain to be allowed to evolve during the iteration process. 

1 0 Figure 6 shows a flowchart illustrating the iterative process of the example embodiment, 

as follows: 

Step 500, applying a normalisation or averaging process to a giv^ FBG amplitude 
function; 

Step 502, solve a direct scattering problem based on the modified grating design 
15 functioI^ 

Step 504, modify the resulting spectral response; 

Stq> 506, solve an inverse scatt^ing problem based on the modified spectral response; 

Step 508, determine whether a desired quahty has been achieved, and 

repeat the steps (retum loop 510) if not, or stop the iteration as indicated by arrow 512. 

20 In the following, a description will be given of how to obtain a suitable initial (pre- 

iteration) grating design (compare Figure 1), in example embodiments. 

In the example embodiment, a multi-<)hannel grating function is constructed by solving 
the standard inverse scattering problem for multi-channel wavelength-shifted spectral 
charactaistics. Spectral response functions Hr (X) for partial single gratings are being de- 
25 phased with respect to each oth^. In other words, the inverse scattering problem mscy be 
presented as: 
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with nonzero relative phases in the example embodiment It is to be noted that 
spectral responses are being de-phased, not partial seeding gratings tiiemselves. After solving 
the inverse scattering problem for Hr (X), the multi-channel grating design function is 
obtained and may be presented in a form: 

5 q{z) = ic,j2^'^^^*^*^^ s K e'^""^^^ (2) 

where we explicitly retain, for illustrative purposes, the single channel grating design 
function The remaining fectors in the expression (2) represent a "sampling^' 

function, which is qjeriodic. It is thus not a sampling, i.e. periodic, function anymore, but 
deviates from a sampling function. 
10 It will be q^eciated by a person skilled in the art that in order to initially (prior to our 

iterative sch^e) pre-optimise the grating design for minimxmn refiractive ind» change An^ as 
a function of the number of channels N, suitable numerical or analytical methods can be 
applied. 

For relatively small number of channels N one can numerically scan through possible 
15 combinations of relative phases i9, (solving an inverse scattering problem for each particular 
combination of dephasing angles Si) and selecting the combmation which is optimal according 
to some specific selection criterion (e.g. selecting the combination which minimises maximum 
required refractive index change). 

For N»l location of the optimal set i9, is diflScult Even rough direct scanning 
20 tihrough all possible sets of angles (followed by efficient numerical miniminn search routines) 
quickly becomes numerically inefScient. 

In the example embodiment, the approximate equivaloice between partial spectra 
dephasing angles &i and^artial grating relative phases ^ is utilised. Indeed, for weak gratings 
the first order Bom approximation holds: 

25 --q(z/2) = rr(J3)exp{--i^)d/3, (3) 

2 ^ 

where 9(2) is a grating design function and r(J3) is a complex reflection coeffideni 
The Fourier transform (3) is a linear operation with a major property 
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15 



20 



pf^j j(0 ^ jfl)^ ^ p /I) ^ ^^p Thus, in this approximation, dephasing of partial 
gratings is equivalent to dephasing of partial spectral channel responses. Formally the last 
statement does not hold beyond weak grating limit. However, in practice, it is still 
approximately correct and the optimal set of angles (for dephasing of partial gratings) may be 
used as a very good approximation for the optimal set of partial spectral channel angles i9/. The 
only important change related to flie "sampling" fimction between dephasing of partial channels 
and dephasing of partial gratings approaches is that it becomes slightly 2?)eriodic, as mentioned 
above. 

Therefore, the phase shift values &i may be taken &om a partial gratings de-phasing 
grating design method In the example embodiment, a sampling function which periodically 
modulates the amplitude of a given single-channel grating (seeding grating) is utilised. In 
addition to tiie periodic modulation of the amplitude of the seeding grating, different relevant 
phases ^ for each of the wavelength-seeding gratings are introduced Accordingly, the 
resulting design fitnction in the exanaple embodimrat may be expressed as: 



(4) 

where the additional phase of the grating yf^yf{z) and the sanipling ampHtude Q=Q(z) 
are given by: 

\HI2 , 



g2(2) = 45]cos(a; -a^jcos(n,AAz/2+p;)cos(n^Afe/2 + p^), and 
^l^l^sina/ cos(n,AAz/2 + >ff^) 



yf{z) = tan"* 



^J^l^cosa; cosififMz/l-^ fi^ 



is even 



or 



(Ar~i)/2 

2^(^) = 4 ]^cosa,cos(n,Afe/2+P;)+ 
/•I 

iN-l)/2 

4 ^cos(a; -a^)cos(7i;Afe/2+p;)cos(7i^Afe/2 + p^)+l, 



and 
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= tan 



-J 



^j^''^"sma, cos(n,Afe/2+p,) 
2^,"""cosa, cos(n,Afe/2+p, +1) 



^ is odd. 



where ^2l-N-l and «^ =2p-N-l. 

Jn the above expressions for and ^z) we use notations 

as($), +^;^+,_,)/2,A -^y^+w)/2andset $}(;,+,)/2 =0 for odd number of channels. 

5 hi Older to reduce Ihe maximum value of Q along the grating structure, while trying to 

avoid touching the zoo level along the grating stracture, the minimization strategy used in the 
example onbodiment is minimizing its maximum deviations of Q(z) along z from the 

theoretical limit level of 7^. Mathematically this may be formulated as finding 

a.(.;ar,A->) forwh^ch ^^^^^^[^3x,{e(.;a.A)}-nun,tea,.A)}l 

10 This approach is iiiq)lemeDted in the cKamplc embodiment by using a simulated 

annealing algorithm - a Monte Carlo approach for minimizing of multi-variable functions. 

The phase shift values (a/''"^^/*^'^) thus determined are thai utilised as phase shift 

values Si for the inverse scattering problem as part of the example embodiment (see equation 
(1)). The obtained dephased ^-channel design may be used as an input for the iterative 
1 5 optimisation procedure anbodying the present invention. 

In the prefen^ embodiment, the strategy of difference minimisation is then extended 
further. The JV-central-channel-only minimisation set of dephasing angles ^ (1 <i<N)is taken 
and used to constmct the initial sampling function S(z) = 0(z) e on a period 2n/M, Then the 
nontrivial Q(z) dependence is replaced by the asyinptotic value 1^^ and the resulting 

20 intermediate sampling function (with only phase sampling preset) is decomposed in a Fourier 
series to jSnd new Fourier coefficients a/ and a new set of ttie corresponding dephasing angles 
s Arg(fli). This procedure leads to a change in partial higji-amplitude grating coefficients. In 
addition, n extra small higher order harmonic terms (corresponding to additional small 
amplitude partial gratings) arise. Next we replace coefficients ai for N high-amplitude partial 

25 gratings with Ce'*^"^' , where a constant C is found from (he normalisation condition 
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tM-NC' (4) 



and all coefficients ai outside the band / e [1 - n / 2, iV^ + n / 2] being set to zero. These modified 
N+n amplitudes are used to construct a new "sanq)ling" function via inverse Fourier 
transformation. 

5 The procedure is repeated until a selected quality criterion has been reached. In the 

example embodiment, the itaation is stopped when a variation in the mayimiiTri "sanq)ling*' 
function amplitude between iterations falls below 10*^. 

The extended dephasing procedure for an infinite number of side channels n translates 
all nontrivial amplitude modulation of the sampling function mto its phase t//(z). As a price for 
10 doing that, about ION additional partial gratings (higher order harmonics with \ai\> 0.001) 
appear ia the spectral characteristics, which means that phase dependence ^z) has a very fine 
structure. However, for practical qjplications the scale of the fine structure should always be 
larger than the size of the laser beam used for grating writing. An estimate for the vspper limit 
for the maximum number iV+n of nonzero partial gratings can be given as 

15 N+n=: — ^ , (5) 

2zoA^«o 

where Ao is the central wavelength [in microns], zo is flie minimal laser beam size [in 
microns], M is the neighbouring channel spacing [in nanometres], and no is the FBG average 
refiractive index. Estimate (5) is essentially based on the Nyquist theorem requiring at least two 
sampling points to define the cosine wave cycle. For the purpose of this desoiption, the xxppcr 
20 limit for the maxiTnum number N+n given by the resolution limitation in the fabrication of the 
grating [estmiate (5)] will be referred to. as un-limited in the context of the preferred 
embodiment of the present invention, as it relates to a given experimental limitation. 

In the example embodimrat, the number of additional small anq)litude partial gratings 
has been (further) limited. It was found that selecting n=2A^+2 channels gives good results in 
25 preferred embodiments of the invention. 

Accordingly, in the example embodiment the initial (pre-iteration) grating design 
function (compare Figure 1) is arrived at through a multi-chaimel grating design approach based 
on partial smgle grating dephasing, which may be further modified by an extended grating 
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dq}hasmg procedure allowing for a limited number of side channels in the spectral 
characteristics. The set of angles and amplitudes obtained by the grating dephasing procedure is 
xised for spectral dephasing of partial channels to obtain a pre-iteration design. After this pre- 
iteration design is obtained, the iterative procedure embodying the present invention is applied. 

5 The implementation of the multi-channel grating design of the preferred embodiment in 

a gratmg structure reqmres grating writing apparatus with high spatial resolution to be utilised. 
Therefore, in a grating writing apparatus relying on photo induced refractive index changes, the 
aqpparatus preferably comprises a beam focusing means to reduce the size of tiie beam in the 
core of the photosensitive waveguide. 

10 Figure 7 shows an example experimental set up SO for writing a multi-channel grating 

52 into an optical fibre 54. The experimental set up 50 comprises an interferometer 56 which 
includes a first acousto-optic modulation 58 being operated under an acousto-optic wave of a 
first frequency Qu as indicated by arrow 14. An incoming light beam 60 is incident on the first 
acousto-optic modulator 58 under a first order Bragg angle. The operating conditions of tiie 

15 acousto-optic modulator 58 are chosen such tiiat the modulator 58 is under driven, whereby 
approximately 50% of the incoming beam 60 is diffracted into a first order beam 62, and 50% 
passing through the acousto-optic modulator 58 as un-dififracted beam 64. Hie un-diffracted 
beam 64 is incident on a second acousto-optic modulator 66 of the interferometer 56 under a 
first order Bragg angle, whereas the beam 62 is not. Accordingly, the beam 62 passing througjh 

20 die second acousto-optic modulator 66 without any significant loss. 

The second acousto-optic modulator 66 is operated under an acousto-optic wave of a 
frequency Qi, which propagates in a direction opposed the direction of the acousto-optic wave 
in the first modulator 58 as indicated by arrow 68. After the second acousto-optic modulator 66 
the first order diffiracted beam 70 and the beam 62 are frequency shifted in the same direction 
25 (e.g. higher frequency), but by dififerrat amounts i.e. Cl\ v 

The beams 62, 70 are tiien brou^t to interference utilising an optical lens 72, and tiie 
resulting interference pattern (at numeral 74) induces refractive index changes in the 
photosensitive optical fibre 54, whereby a refractive index profile, i.e. grating structure 52, is 
created in the optical fibre 54. 
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In Figure 7, the optical fibre 54 is translated along the interferometer at a speed v, as 
indicated by arrow 74. 

It will be q)preciated by a person skilled in the art that the experimental set SO shown 
in Figure 7 can be utilised to write a multi-channel grating structure of a multi-channel grating 
S design embodying the present invention througji suitable control of the first and second acousto- 
optic modulators 58, 66, in conjunction with a suitable control of the speed v at which the 
optical fibre 54 is translated along the interferometer 56 at any particular time. The high spatial 
resolution reqmred to implement the multi-channel design of the prefenred embodiment is 
achieved in the set vp shown in Figure 7 by utilising optical lens 72, with the practical limit of 
10 the beam size in the focal plane preferably being of the order of the waveguide core size. 

It will be qspreciated by a person skilled in the art that nummus variations and/or 
modifications may be made to the present invmtion as shown in the specific embodiments 
witibiout departing fit>m the spirit of scope of the invration as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be iUustrative and not restrictive. 

15 For example, it will be appreciated that the initial grating design fimction utilised as the 

starting point of an iniprovement method embodying the present invention may be derived 
utilising different methods as the ones described for the example embodiment. Alternative 
methods can be found e.g. in Patent Co-optation Treaty (PCT) patent ^plication no. 
PCT/AU02/00160 filed on 15 February 2002, entitled 'Multi-channel grating design" assigned 

20 to the present ^licant, PCT ^fication no. PCT/AU02/00309 filed on 15 March 2002, 
^titled ''Grating design" assigned to the present applicant, and/or Australian provisional patent 
plication no. PS1131 filed on 15 March 2002, ^titled "Improved multi-channel grating 
design" in the name of the present applicant. 

Furthermore, it will be ^predated that the present invention is not limited to the type of 
25 multi-channel gratings of the example embodiments described. For example, the present 
invention can be equally applied to multi-channel gratings consisting of non-identical spectral . 
chaxmels, and/or gratings consisting of groins of spectral channels. 

Furthermore, it will be appreciated that otibier numerical or mathematical approaches can 
be taken to implement the iterative procedure embodying the present invention, e.g. comprises 
30 solving complex integral equations. 
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Furthennore, multi-channel gratings can be created on the basis of the mnlti-channel 
grating design of the present invention using various known grating creation techniques, 
including one or more of the group of photo-induced refractive index variation in photo 
sensitive waveguide materials, etching techniques including etching techniques utilising a 
S phasemask, and epitaxial techniques. 

Also, it will be £^preciated that methods embodying the present invention may be 
complimented by known optimisation techniques such as generic optimisation algorithms. 

In the claims that follow and in the summary of the invention, excq>t where the context 
requires otherwise due to express language or necessary implication the word ^'comprising" is 
10 used in the sense of "including', i.e. the features specified may be associated witii further 
features in various embodimmts of the inv^otion. 



